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Ultra-fast x-ray tomography for multi-phase 
flow interface dynamic studies 



The present paper describes the concept of* fast wanning X~ 
ray tomograpJi, the hardwire development, and tnewutremetu 
results of gas-liquid two-phase flow in u vertical pipe. The de- 
vice uses tS pulsed X-ray sources activated in a sueeesiive or- 
der. In this way, a complete set of JS independent projections 
of the object is obtained within 3$ ws, L e. the measuring rate 
is about 250 frames per second. Finally, to evaluate the mea- 
surement eapabillxy of the fust X-my C% a wire-mesh sensor 
was installed hi the flow hynp and both systems were operated 
for the same two-pfutse flow sirmdhmeousty. Comparison of 
rite time series of the eross section averaged void fraction from 
both systems xftowed sufficient agreement for shtti pow at large 
void fractions, white the fust CT underestimated the void frac- 
tion of bubbly flaw especially in law void fraction rttnge For 
the wire-mesh sensor, curved deformation of slug bubble in- 
terface was found, further luudmoe improvement is in pro- 
gress to achieve tetter resolution with the fast X-ray CT scan- 
ner. 

Vltra-schnelle Hontgentomographie ftir das Stadium derOy- 
nam Ik von Phaxengrtnzfladten in Mehrphasenstrumungcn. 

Der He.itrag beschreibt das Konzept eines schnelfen Jiontgtnto- 
mographen, die llardwarcentwivklnng und Messerge.bnHse am 
einer Zyveiphusenstrarmotg tines Flussigkeiis-Gas-Geuvschs in 
einer vertikalen RtiksieiUtfifr Das Oeriii verfiigt iiber 18 ge- 
puhte ROmgenqurlten, die m schneller folge aktiviert v*erdr.ri. 
Auf diese Weise wird em vollstandiger Soft aus 18 Projektio* 
rten des Objekts inrwrhalh von 38 mt gewonnen. Die Messfre- 
(jitenz betragt demznfolge «i. 250 BiUer pro Sekunde. Urn 
das System zu test en. n-urtle. ein Gittersensor in die Stro- 
mungsschieife eingehaut HeMe Cerate warden gemeinsatn 
zur Auf/uditne der gleichen Zweiphasenstrbtmmg genuizt 
Der Vergieich des Ztitverlaufs des ttuersctmittsgemittetten GtiS- 
gehutts zeigt ewe gate Oberebtsiimmung z* ! isdten beiden Ver> 
fnhren im Gebiei der f*fropfen#trvmung. Bei einer Masenstrti- 
mung komtnt es zu enter Unterbewemtng des Gasgehate 
durch den Hontgentontographen. In Bezug auf den Ginersen- 
sor veurtle eine Deformation der Grvflblasen gefunderu Zu- 
kOnfiige Hardware-Verbesserungen zur Erhohitng der Aufltf- 
sung des HUmgemomographen siml geplam. 



1 Introduction 

l.Tse of X-ray CT (Computer 'tomography) system has been 
extended from medical fields to a variety of industrial appli- 
cations [1, 21- Among them is the visualization of muh {-di- 
mensional void distribution of two-phase ilovt in lube?. Mitsu- 
lake et al. [3] and Mornoka ci al. [4| used medical X-ray CT 
scanners modified for void fraction measurement within rod 
bundles. The results vvere comparable with (heir predictions 
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based on three-Quid model Recently Tiseanu et al. applied 
cone bc«rn lomograpliy to bubbly flow in a vertical pipe and 
obtained Sharper void peaks near the wall than previously re- 
ported [5J. One beam iwnc^raplty is a 3D visualization lech- 
nique and allowed them to collect raulli-sllce data in vertical 
direction within the same amount of scanning time as 213 sy*- 
lenvt Gamma-my is another form of geiierating hi^h energv^ 
photons. Kumar et aL bvalc a tomography yystem based on 
0-137 gamma ray source for void distribution measurcmeat 
of a btibbte coluina |6J. Although gumma ray sources eliuii- 
naie the beam-hardening effect, lew photon flux results in 
long scanning lime. However, these systems arc not applic- 
able to dynamic interface detection of muhi-pliase flow be- 
came oi low Time resolution compared to transition time of 
interlaces. To visualize dynamic interfaces in multi-phase 
flo«\ time -resolved X-ray CT scanners have been developed 
by Mori and Akai |7], Hon c t al. and Misawa et al. |9| 
(Fig. 1). Scanning time of Ihese s>-stems is reduced by an or- 
der of two compared to existing ones. 7n pursuit of time reso- 
lution, however, spatial resolution is sacrificed due lo less 
numbers of detector elements a>id views. The present paper 
describes the concept of fasi scanning capability, hardware de- 
velopment, and measurement results of gas-liquid two-phase 
How in a vertical pipe- Finally, to evaluate the measurement 
capability of the fast X-ray CT the wire-mesh sensor [10] 
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was installed in the flw loop and both systems were operared 
for Ihe same two-phase flow simultaneously ~I*he results 
showed that the time series of cross section averaged void 
fraction agreed to each other lor sing flow at large void trac- 
tion, but lhat the wire-moth sensor was suitable for bubbly 
flow al low void fraction, l-rom the comparison oC the time 
series data, coerced deformation of slug bubble interlace was 
suggested for the wire-mesh sensor. 



2 Selection of parameters based on simulation 

The requirements set for the multi-phase flow measurement 
were 1) spatial resolution < 3.0 mm, 2) time resolution 
<A msec, per scan, ami 3) measurement area diameter 
<5Drnm, Before hardware construction, numerical simula- 
tions were conducted to see whether the above requirements 
were achieved. Cylindrical and spherical plastic phantoms as 
well as their combination were numerically modeled and their 
projection data was calculated based on the X-ray energy, at- 
tenuation coefficient of the materials, distance between the 
X-ray sources and detectors. Among several key parameters, 
the number of X-ray sources and detector pixels are most im- 
portant parameters which determines uniquely the spatial re- 
solution in the reconstructed image* Fig. 2 shows the recon- 
structed images of a cylindrical rod with several vertical 
holes of different diameters. Spatial resolution varies by chan- 
ging the number of sources find the number of detector pixels. 
In fact, current system consists of IS X-ray sources and 60 pix- 
els per fan beam, which tteorcticaliy allows us to observe 
2 mm diameter boles under a situation free of noise and misa- 
lign mem. One can see that an increase in the X-ray sources 
arc mure efficient to improving die resolution than lhat in 
the pixels in this design, However, as the X-ray sources in- 
creases, the time resolution and intensity will suffer because 
relatively large diameter of the X-ray tubes results in long dis- 
tance between X-ray sources and detectors. 



J J System outline 

To complete a single scan on the order ol milliseconds, the 
mechanical scanning motion should be avoided. Instead, the 
18 X-ray sour ces are placed along a circle of radius 120.5 mm 
from the center ol measurement region, and emission of X- 
ray is controlled by electrical switching. Fig. 3 illustrates the 
switching mechanism ol the X-ray tubes, l^leerrons Keep 
emitted from the healed filament of electron guns but al- 
lowed 10 pass through the bias grid aperture only *hen the 
negative bias was removed. After passing through the grid 
aperture, the electrons are accelerated in the 40-160 kV elec- 
tric field and hit die tungsten target embedded in the common 
anode ring. The X-ray fan beam spreads over 24 degrees and 
activates 6t> pixels on the detector. High voltage is continu- 
ously applied on the anode ring from the voltage source out- 
side" the vacuum chamber. All x-ray sources without glass 
tubes are installed in the vacuum chamber. This configuration 
was adopted to reduce the pitch of X-ray focuses for increas- 
ing spatial resolution. To avoid the interaction of X-ray beams 
with the detector. X-ray focus plane is offset from the detec- 
tor plane by 10-20 mm. This allows for photons to transmit 
tluough the object and reach the detector elements on the op- 
posite side of the X-ray source at an elevated angle. From the 
point of image reconstruction, data collection along a path 
elevated from the horizontal plane can he a cause ol artifact. 
However, simulation result showed that the effect is not sig- 
nificant if the elevation angle is less than 15 degree. f* mission 
Of X-rnv and data acquisition is synchronized by a pulse con- 
troller. Following a start trigger, the bias grids are sequen- 
tially disarmed for 100 microseconds and all the detector ele- 
ments are scanned upon each X-ray emission. During the 
next 100 microseconds, the X-ray sources rest and the data is 
transferred to the computer. After the IS'** X-ray emission is 
done, another data acquisition is made without X-rays. This 
serves at a oflset level for the data sequence in a single scan, 
prom a total of 19 data set. thus acquired at different angles, 
a cross suction image is reconstructed. Table 1 lists the main 
design parameters of the fast X-ray CT scanner. With the gi- 
ven parameters, * spatial resolution of 2.6 mm and the scan- 
ning rate of 263 list can be attained. Image reenrtst ruction, 
base on a filtered back projection, is carried out on off-line 
after all the data is stored in a tiff- type format. The littered 
back projection method is used for 2D image reconstruction 
in 101 x 101 pixel format over a 50 mm diameter region of in- 
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Tabic J. Fast X-nry CT speciftcitiiorut 



Number <\f X-ray sources 


IK 


AcccJcr atJon voltage 


160 fcv max. 


Deani current 


IOioA max. 


$c»n lime per slice 


3& msec 


Scan rate 


263 slices 


Ream jhiclcnes* 


5 mm 


Number of detectors 


256 


Detector device 


CdTe 


No. at dctcctorvYim beam 


60 


Maximum object diameter 


50 mm 


Designed spatial resolution 


2,6 mm 


Fan beam angle 


24 degrees 


X-ray source ring diameter 


241 mm 


I'Jetecrar ring diameter 


1*3.4 mm 


A»D t evolution 


12 bit 


Integration time 


m*cc min. 


Image reooitsu-uctioa 


FBP 


Filler Type 


Shepp-Ivsjgarj 



terest. A reference wale of 50 mm was digiti?ed into 1 01 piv- 
els converting ro 0.495 mm in length or 0:245 mm 2 in area 
per pixel. 



3.2 Devekfprmnu of a high effuitney semi-conductor dtUtcivr 

How to collect a sufficient amount of photon* In the detector 
during the short period of pulsed X-ray irradiation is the key 
issue of detector design. Semi-conductor solid detectors con- 
vert incident photon energy to carrier signals directly inside 



the crystal, resulting in higher conversion efficiency than scin- 
tillator detectors by an order (Fig. 4). Cadmium lelluride 
(CdTe) crystal was selected as the detector material because 
of its large atomic number and large band gap. For the X-ray 
energy around 100 W t more than 95 % photons arc captured 
by a 1.0 mm thick crystal. The band gap of CdTe crystal is 
1.44eV, arid this allows us to use the device without cooling. 
On the surfaces of a substrate of CdTe monocrystal (size: 
5 x 20 x 1 mm), anode and cathode pixel patterns arc printed 
by pholohthography technique. The effective pixel size is 
0.3 x0,8 mm with « pitch of 1.425 ram. A set of 32 pisels is al- 
located in a single detector module. Signal processing of this 
charge-sensitive device is made by a 32 channel CMOS IC. 
The chip is placed close to the sensing crystal to reduce the 
noise effect. The developed detector consists of 8 modules 
having 32 pixels each and completes scanning over 256 pixels 
in 200 micro seconds. Alter a successful radiography tests* 
the detector modules were assembled into the fast CT scan- 
ner, 



4 Dynamic interface measurement 

4. 1 Examination of cUifa qtutttiy 

Prior to fluid experiment, measurement capability of the 
cross- seel ion and the interface was evaluated by scanning 
plastic phantoms with know dimensions. In I he first scries 
of experiment, spheres and cylinders of different size* 
moved through the measuring plane and their volumes and 
surface areas were compared with true values by setting a 
threshold value based on intensity histogram. T!*s results in- 
dicate that the volumes and interfaces for velocities less 
than t.Onvs can be reproduced within 20 % of mcasurcrneni 
accuracy. 

Nexl. a slug unit phantom was made to simulate more rea- 
listic situation. Hie phantom has a semi-spherical dome, a cy- 
lindrical part of 2 mm wall thickness, and cylindrical and 
spherical voids of mm in diameter. The reconstructed 
images of several cross sections for the moving velocity of 
0.3 nv's arc shown in fig, 5. When the dome passes through 
the measuring plane, its bound tiry looks diffused because the 
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Fi^. * Hr.cinnMrm:ied image* of a aluv shape phantom 

cro^> section shape was constantly changing. » n the axial direc- 
tion while being scanned. On conirary, the cylindrical part 
was reproduced clearly since the cross section shape rental ncd 
the same during the scanning period. Cylindrical and spheri- 
cal cavities of 3-6 mm in diameter were also reproduced suc- 
cessfully. 

fig. 6 shows the average intensity profiles of reconstructed 
iraa&es for water a lone in a 42 mm diameter pipe and for bub- 
bles in the water. Air was introduced through a 4 mm diameter 
nozzle placed at the bottom ol the pipe The average intensity 
fluctuation for water is considered to be measurement noise. 
When a signal (ally below this level, identification of individual 
bubbles becomes difficult. However, from the fact thai Ihe 
noise signal fluctuates around level aero, the average void frac- 
tion over the cross section can be measured lrom level zero 
with an error on the order of noise fluctuation. 

4. 2 Air- water two-phase jlo w 

With the knowledge of signal characteristics in measurement, 
the superficial gas velocity was increased from 0.013 to 
Q.091 m/s. For visualization of interlace, the attenuation coef - 
fjciem matrix ts mapped to in the 256 grayscale level. A 
threshold level was chosen by referring to histograms for each 
condition. The figures shown in Fig. ? represent the time ser- 
ies of spatial void distribution at the measuring plane. Recoti* 
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structed tmuges in W X101 pixel formal are stacked in order 
with a slice thickness of one pixel. Tn reality, twwever, the lo- 
cal interface velocity are dunging in time. Ihe best approxi- 
mation available is to use the average bubble velocity after 
calculating the average void fraction. The overall distribution 
patterns are sufficiently reserved regardless o£ instantaneous 
fluctuation of inter faces 

In addition to slug flow, bubbly flow was created by intro- 
ducing air through a porous metal with its pore size of 100 mi- 
cro meter. Although the individual bubble shape is difficult to 
recognize, the concentration of void in the cross section is vi- 
sualized. The image shows the concentration of bubbles near 
the pipe wall for concurrent How under relatively low gas flow 
rate*. As the gas flow rate increases, bubble concentration 
starts to occur in the central region of the pipe 

4.3 Evaluation of dam 

The measured void distribution in the cross section contains 
information on the void fraction and the interlace area con- 



Fig. 6. Average intensify profile for wler tdonr ami txthbirs in n>aiet 
f/jr a (1091 m/> ) 
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centration, both of which arc essential quantities ia multi-di- 
inensianal two-phase flow analysis. After image reconstruc- 
tion, a slice disk including 2D bubble distribution is obtained 
S). Assuming that the void distribution is steady in time 
And the slice thickness is sufficiently thin, the void fraction is 
the ratio of the total bubble cross section area Sum(A^i) to 
the total flow area 

a«I(A*yAg II) 

Oi> the other hand, the interfneia! area concentration is by de- 
finition the amount of interface in a unit volume. The amount 
of interface is calculated as ihe product of the bubble bound- 
ary and Ihc thickness of the disk. Thus, 

Since the total flow area remain unchanged, both quantities 
depend on the cross section and perimeter of the bubbles in 
the slice disk. 



15mm 




f*i>, 9. Wire-mcih \<n.wr and Ihe (tai X-my CT measurement partes 

1052 slices. The superficial gas velocity ranged from 0.02 to 
0,*10 m's. while the liquid velocity varied from 0.0-0.69 m/s, 
lioth bubbly and slug flow patterns were tested under the 
same flow conditions. 



5 Comparison with the wire-mesh aeo&or 

Advantages of ihc last X-ray CT are its non-invasivcness 
and applicability to all types of mu hi -phase How including 
solids. Wire- mesh tensor is especially useful when an on- 
line monitoring is necessary and the flow velocity is rela- 
tively large. However, it is difficult to calibrate the mea- 
sured data each alone since the absolute values are often 
unknown in actual two-phase flow. When both systems are 
applied to the same two-phase flow, it is curious to see 
how much the instantaneous values differ from each other. 
Both methods were applied to Che same air- water two- 
phase How to validate their measurements as well as to 
clarify the difference, 

i.l Experimental scrap 

lite 16x 16 wire-mesh sensor and the fast X-ray CT scanner 
were applied simultaneously to the two-phase flow* ia a 
42 mm diameter vertical pipe. A part of the apparatus U 
shown in Fig. 9. Since the fast X-ray CT is non-invasive, its 
measurement plane was set at 15 mm upstream of that of 
ihe wire -mesh sensor. This interval w as the minimum to 
avoid interference of each other's measurement. The mea- 
suring plane was located 1.15* m from tlie air-water mixing 
chamber. The starting trigger was shared by the both sys- 
tems and the darn for 4.0 s. were recorded under various 
flow conditions. Sampling frequency for the wire mesh sen- 
sor and the fast X-ray CT were set at 1053 H*/. and 263 Hz. 
respectively, the total cross section images for the wire- 
mesh sensor was 4211, while the fast X-ray CT recorded 



5,2 Jtcsufo of voitJ fmahm measurement 

With the wire- mesh sensor and the fast X-ray CT, large slug 
bubbles were similarly reproduced, except that tlie interface 
shape appears to be elongated for the wire-mesh measure- 
ment. This is caused either by the interface repulsion against 
the sensor wires or forsnaiion of liquid bridges between the 
sensor wires. Hi^h -speed video observation during the wire- 
mesh measurement also indicated such cases may occur. On 
the other hand, the fast CT measurement provided relatively 
small values for the void fraction component in the liquid 
hold up. This difference comes from Ihe sensitivity of the 
measurement system toward small bubbles. As the invasive 
method, the wire-mesh sensor directly senses the existence 
of small bubbles while the X-ray contrast becomes smaller 
as the bubble si&e is reduced In processing of the fast X- 
ray CT data, a constant threshold level was used to make 
binary images for void fraction count. Use of variable 
thresholds depending ou the bubble size will improve the 
measurement 

Tune series data of the cross section as r e raged void fraction 
were compared for the both measurements. The wire -mesh 
data was re-sampled at 263 Hx to match the sampling fre- 
quency of the last CT. Then the average time lag was ob- 
tained by calculating the cross correlation of the profiles, 
fig. 10 show? the void fraction time series after re-sampled 
and adjusted by the time lag. Again the slug bubble profiles 
are similar to each other, except the front -edge deformation 
as seen in the reconstructed image. Histograms Of the void 
fraction values for both measurements indicates their differ- 
ence clearly; Although the overall shape resembles, the peak 
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Fig ft. Cvmparis&x of cross Section average void frarttnn mtramrtd by 
the wire-mesh s&tt&rimtl the fast X-ray CT 

of the vvirc-mesh measurement near void fraction 0.5 may be 
caused by false pealr due to coerced deformation of slug bub- 
bles. On the other hand, the large peak near void fraction 
zero in the fast X-ray CT measurement indicates its Inabilhy 
to detect small bubbles which the wire-mesh sensor success- 
fully captured. Regarding bubbly How, both measurements 
showed similar values for flic cases of large gas flow rales. II 
is difficult 10 resolve single hubbies in this flow regime. When 
l he number of hubhle* in the liquid becomes less al low gas 
flow rates, the chance of capturing no bubble region will in- 
crease for the fast X-ray CT, resulting in small values of aver- 
age void fraction, In ihis sense, the wire- mesh sensor with 
high sampling licqueney is suitable for bubbly flow of low 
void fraction, Fig- 1 1 shows the comparison of the average 
void fraction obtained from the wire-mesh sensor and the fast 
X-ray OF for bubbly and slug flows. A tendency is observed 
thai the fast X-ray CT under estimates the void fraction espe- 
cially in bubbly liosv regime. To improve the measurement 
quality for the fast X-ray CP, the following measures arc un- 
dertaken: 1) adaptive selection of threshold for bmarizaiion. 
2} development of detector with largc-siax pixels, and 3) de- 
velopment of X-ray sources of uniform intensity. 



6 Conclusion 

Concept and operation of a fust X-ray CP scanner were de- 
scribed in detail. A highly sensitive X-ray semt-coiuJuctor dc- 
teelor was newly developed and installed to increase the de- 
tection capability of the system. Alter the temporal and 
spatial resolutions were examined by both simulations and 
the phantom tests, the system was applied to air- water two- 
phase flow. Detailed dynamic interlace was successfully visua- 
lized by the system. 

Tb verify the interface measurements, the wire-mesh sensor 
was installed close to the measurement plane of the fast X-ray 
CT. The same two-phase flows were simultaneously measured 
by the both systems and their results were compared to each 
other. Comparison of the time series of the cross section aver- 
aged void fraction from both systems showed sufficient agree- 
ment for slug flow at large void fractions, while the fast CT 
underestimated the void fraction of bubbly flow especially in 
low void fraction range. Further hardware improvement is in 
progress lo achieve better resolution with the fast X-ray CT 
scanner. 
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Nomenclature 



Sign 


unit 


Denomination 


ct 




void fraction 


a. 


i;m 


Interfacial area concent ration 


A„ 


ror 


cross section area of a pipe 




m~ 


cross section area of a babble 


Ui 


m 


Interlace length of a bubble 


d* 


m 


slice thickness 


x, y, si 


in 


co-ordinates 


h 


m/s 


gas superfieial velocity 


h 


m/s 


liquid superficial velocity 
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